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We treat the problem of diffusion of solute atoms around screw dislocations. In particular, we 
express and solve the diffusion equation, in radial symmetry, in an elastic field of a screw dislocation 
0^ | subject to the flux conservation boundary condition at the interface of a new phase. We consider an 

incoherent second-phase precipitate growing under the action of the stress field of a screw dislocation. 
The second-phase growth rate as a function of the supersaturation and a strain energy parameter 
is evaluated in spatial dimensions d = 2 and d = 3. Our calculations show that an increase in the 
, amplitude of dislocation force, e.g. the magnitude of the Burgers vector, enhances the second-phase 

■ growth in an alloy. Moreover, a relationship linking the supersaturation to the precipitate size in 

Pj^ ' the presence of the elastic field of dislocation is calculated. 
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I. INTRODUCTION 



Dislocations can alter different stages of the precipitation process in crystalline solids, which consists of nucleation, 
growth and coarsening p| Q • Distortion of the lattice in proximity of a dislocation can enhance nucleation in several 
ways 0, H| • The main effect is the reduction in the volume strain energy associated with the phase transformation. 
Nucleation on dislocations can also be helped by solute segregation which raises the local concentration of the solute 
in the vicinity of a dislocation, caused by migration of solutes toward the dislocation, the Cottrell atmosphere effect. 
When the Cottrell atmosphere becomes supersaturated, nucleation of a new phase may occur followed by growth 
of nucleus. Moreover, dislocation can aid the growth of an embryo beyond its critical size by providing a diffusion 
passage with a lower activation energy. 

Precipitation of second-phase along dislocation lines has been observed in a number of alloys 0, For example, 
in Al-Zn-Mg alloys, dislocations not only induce and enhance nucleation and growth of the coherent second-phase 
MgZn 2 precipitates, but also produce a spatial precipitate size gradient around them @,H, Hj]. Cahn [l(| provided the 
first quantitative model for nucleation of second-phase on dislocations in solids. In Cahn's model, it is assumed that 
a cross-section of the nucleus is circular, which is strictly valid for a screw dislocation [l|. Also, it is posited that the 
nucleus is incoherent with the matrix so that a constant interfacial energy can be allotted to the boundary between the 
' new phase and the matrix. An incoherent particle interface with the matrix has a different atomic configuration than 
that of the phases. The matrix is an isotropic elastic material and the formation of the precipitate releases the elastic 
t-H ' energy initially stored in its volume. Moreover, the matrix energy is assumed to remain constant by precipitation. In 
this model, besides the usual volume and surface energy terms in the expression for the total free energy of formation 
of a nucleus of a given size, there is a term representing the strain energy of the dislocation in the region currently 
Q\ [ occupied by the new phase. Cahn's model predicts that both a larger Burgers vector and a more negative chemical free 
■ energy change between the precipitate and the matrix induce higher nucleation rates, in agreement with experiment 

> : Iff. 

•rH , Segregation phenomenon around dislocations, i.e. the Cottrell atmosphere effect, has been observed among others 
in Fc-Al alloys doped with boron atomsQjJ and in silicon containing arsenic impurities [l^] , in qualitative agreement 
with Cottrell and Bilby's predictions [13 |. Cottrell and Bilby considered segregation of impurities to straight-edge 
dislocations with the Coulomb- like interaction potential of the form </> = Asm9/r, where A contains the elasticity 
constants and the Burgers vector, and (r, 6) are the polar coordinates. Cottrell and Bilby ignored the flow due 
concentration gradients and solved the simplified diffusion equation in the presence of the aforementioned potential 
field. The model predicts that the total number of impurity atoms removed from solution to the dislocation increases 
with time t according to N(t) ~ i 2 / 3 , which is good agreement with the early stages of segregation of impurities 



to dislocations, e.g. in iron containing carbon and nitrogen |14(. A critical review of the Bilby-Cottrcll model, its 
shortcomings and its improvements are given in pT 
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The object of our present study is the diffusion-controlled growth of a new phase, i.e., a post nucleation process in 
the presence of dislocation held rather than the segregation effect. As in Cahn's nucleation model [ld| , we consider 
an incoherent second-phase precipitate growing under the action of a screw dislocation field. This entails that the 
stress field due to dislocation is pure shear. The equations used for diffusion-controlled growth are radially symmetric. 
These equations for second-phase in a solid or from a supercooled liquid have been, in the absence of an external 
field, solved by Frank [l6[ and discussed by Carslaw and Jaeger [l7| ■ The exact analytical solutions of the equations 
and their various approximations thereof have been systematized and evaluated by Aaron et al. [l8| , which included 
the relations for growth of planar precipitates. Applications of these solutions to materials can be found in many 
publications, e.g. more recent papers on growth of quasi-crystalline phase in Zr-base metallic glasses [l9j and growth 
of Laves phase in Zircaloy [2(| ■ We should also mention another theoretical approach to the problem of nucleation 
and growth of an incoherent second-phase particle in the presence of dislocation field [2l|. Sundar and Hoyt [2l| 
introduced the dislocation field, as in Cahn (To| . in the nucleation part of the model, while for the growth part the 
steady-state solution of the concentration field (Laplace equation) for elliptical particles was utilized. 

The organization of this paper as follows. The formulation of the problem, the governing equations and the formal 
solutions are given in section [ill Solutions of specific cases are presented in section IIII| where the supersaturation 
as a function of the growth coefficient is evaluated as well as the spatial variation of the concentration field in the 
presence of dislocation. In section HVl besides a brief discourse on the issue of interaction between point defects and 
dislocations, we calculate the size-dependence of the concentration at the curved precipitate/matrix for the problem 
under consideration. We have carried out our calculations in space dimensions d — 2 and d — 3. Some mathematical 
analyses for d = 3 are relegated to appendix [Al 



II. FORMULATION AND GENERAL SOLUTIONS 



We consider the problem of growth of the new phase, with radial symmetry (radius r), governed by the diffusion 
of a single entity, u = u(r,t), which is a function of space and time (r, t). u can be either matter (solvent or solute) 
or heat (the latent heat of formation of new phase). The diffusion in the presence of an external field obeys the 
Smoluchowski equation [22[ of the form 

J = D(Vu-PFu), (2) 

where D is the diffusivity, f3 = 1/fcsT, fcs the Boltzmann constant, T the temperature, and F is an external field of 
force. The force can be local (e.g., stresses due to dislocation cores in crystalline solids) or caused externally by an 
applied field (e.g., electric field acting on charged particles). If the acting force is conservative, it can be obtained from 
a potential cj) through F = — V(/>. The considered geometric condition applies to the case of second-phase particles 
growing in a solid solution under phase transformation [20j or droplets growing either from vapour or from a second 
liquid [161 ]. A steady state is reached when J = const. = 0, resulting in u = uo exp (—/?</>). 

Here, we suppose that the diffusion field is along the core of dislocation line and that a cross-section of the precipitate 
(nucleus), perpendicular to the dislocation, is circular, i.e., the precipitate surrounds the dislocation. Furthermore, we 
treat the matrix and solution as linear elastic isotropic media. The elastic potential energy of a stationary dislocation 
of length I is given by [H Hi] 

r 

4> = A In — , for r > ro (3) 

ro 

where A — Gb 2 l/An for screw dislocation, G is the elastic shear modulus of the crystal, b the magnitude of the Burgers 
vector, v Poisson's ratio, and ro is the usual effective core radius. Also, we assume that the dislocation's elastic energy 
is relaxed within the volume occupied by the precipitate and that the precipitate is incoherent with the matrix. Hence 
the interaction energy between the elastic field of the screw dislocation and the elastic field of the solute is zero. In 
the case of an edge dislocation and coherent precipitate/matrix interface, this interaction is non-negligible. 

We study the effect of the potential field ([3]) on diffusing atoms in solid solution using the Smoluchowski equation 
(fTj). The governing equation in spherical symmetry, in d spatial dimension, with B = (3 A, is 

1 du d 2 u . , . 1 du , , . „ u 

n^7 = TT + (d-l + B -— + d-2 B-. 4 
D at ar z r or H 



Making a usual change of variable to the dimensionless reduced radius s — r/y/Dt, the partial differential equation 
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(U|) is reduced to an ordinary differential equation of the form 

d 2 u fs d—l + B\du u 

d^ + (2 + ^^)dI + (d - 2) ^-°' (5) 

with the boundary conditions, zt(oo) = u m , and u(2A) = u s , where u m is the mean (far-field) solute concentration in the 
matrix and u s is the concentration in the matrix at the new-phase/matrix interface determined from thermodynamics 
of new phase, i.e., phase equilibrium and the capillary effect. Moreover, the conservation of flux at the interface radius 
R = 2X\[~Di gives 

K d R d - 1 \3\ r=R = q^, (6) 

where K d = 2n d / 2 /T(d/2), T(x) the usual T-function, V d = 2n d / 2 R d /dT(d/2), and q the amount of the diffusing entity 
ejected at the boundary of the growing phase per unit volume of the latter (new phase) formed. In s-space, equation 
([6]) is written as 

The boundary condition u(2X) = u s and equation |(7J| will provide a relationship between u s and u m through A. 
For d = 2, equation {5} is very much simplified, and we find 

(Bu m + 2q\ 2 )X B e^r(-B/2,s 2 /A) 

u(s) = Um + rWrfflp) — ' (8) 

where we utilized u(oo) = u m and equation l|7|). Here T(a, z) is the incomplete gamma function defined by the integral 
r(a, z) = f°° t a ~ 1 e~ t dt [25]. The yet unknown parameter A is found from relation (jHJ) at u(2A) = u s for a set of input 
parameters u s , u m q, and B, through which the concentration field, equation ([§]). and the growth of second-phase 
(R = 2XyDt) are determined. 

Let us consider the case of d = 3, that is assume that the potential in equation ([3]) is meaningful for a spherically 
symmetric system. In this case, for B ^ 0, the point z = is a regular singularity of equation ([5]), while z = oo is an 
irregular singularity for this equation, see appendix YK\ for further consideration. Nevertheless, for d = 3, the general 
solution of equation ([5]) is expressed in the form 

„ (8) . 2Cl I; I±J> 4) 8 - + 2 * ft - f ; i^,-^.-, (9) 

where ii 7 ' 1 (a; &; z) is Kummer's confluent hypergeomtric function, sometimes denoted by M(a,b,z) (25j. The inte- 
gration constants C\ and Ci in equation ([9]) can be determined by invoking equation ([7]) and also the condition 
u(oo) — u rn , cf. appendix [K\ 

III. COMPUTATIONS 

To study the growth behavior of a second-phase in a solid solution under the action of screw dislocation field, 
we attempt to compute the growth rate constant as a function of the supersaturation parameter k, defined as k = 
(u s — u m )/q u with q u = u p — u s , where u p is the composition of the nucleus [l^|. For d = 2, i.e., a cylindrical 
second-phase platelet, equation ((5J) with u(2X) — u s yields 



k 



2X 2 + Bu m (u p - u s ) 1 



X B e x T(~B/2,X z ). (10) 



2-BX B e^T(-B/2,X 2 ) 

For B = 0, the relations obtained by Frank [l6| are recovered, namely 

u{z) = u m + q u X 2 e x2 E 1 (z 2 /i), (11) 
k = AV 2 £i(A 2 ), (12) 

where Ei(x) is the exponential integral of order one, related to the incomplete gamma function through the identity 
E n (x) =x n - 1 T(l-n,x) [H. 
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From equation (|10p , it is seen that a complete separation of the supersaturation parameter k = (u s — u m ) (u p —u 8 ) 1 
is not possible for B ^ 0. However, for u s << u p (a reasonable proviso) we write 



A 2 + | £ ) A B e A- r 



with e = u s /u p . For B — 1, equations J} 

u(z) 



B/2,X 2 ) + 0(e 2 ), 
and (fl"3|) yield, respectively 

2Ae A2 (u m + 2(7 u A 2 )^ 3/2 (z 2 /4) 



[2-e A2 i?3/ 2 (A 2 )]z 



fc = (^ + -je A S 3/2 (A^) + 0(e^). 



Similarly for B — 2, we have 



(A 2 



4A 2 e A > m + g M A 2 )£ 2 (z 2 /4) 

[1 - e x2 E 2 (X 2 )]z 2 
e)E 2 (X 2 ) + 0(e 2 ). 



(13) 

(14) 
(15) 

(16) 
(17) 



We have plotted the growth coefficient A = R/2\fDt as a function of the supersaturation parameter k in figure Q] 
and the spatial variation of the concentration field in figure [5] for d = 2 and several values of B. The computations 
are performed to 0(e 2 ) with e = 0.01. Figure [T] shows that A is an increasing function of k; and also, as B is raised 
A is elevated. This means that an increase in the amplitude of dislocation force (e.g., the magnitude of the Burgers 
vector) enhances second-phase growth in an alloy. 

Figure [2] displays the reduced concentration versus the reduced radius z = r/yDt for A = 1. The reduced 
concentration is calculated via equation (J5J) . It is seen that for z < 1.6 the concentration is enriched with increase in 
B, whereas for z > 1.6, it is vice versa. So, for A = 1, the crossover z-value is z c w 1.6. Also, as A is reduced, z c is 
decreased. 
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FIG. 1: Growth coefficient A as a function of supersaturation k at various levels of dislocation force amplitude B for a circular 
plate (d = 2) and u s = 0.01m p . 



For d = 3, i.e., a spherical second-phase particle in the absence of dislocation field [B = 0), we find 

-2e- z2 /4 



3„A^ 



u(z) = u m + 2q u \ e 



- V^erfc(z/2) 



k = 2A 2 1 - V^Ae A erfc(A) . 



(18) 
(19) 
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FIG. 2: Reduced concentration field as a function of reduced distance from the surface of the circular plate (d = 2) at various 
levels of dislocation force amplitude B and at A = 1. 



This corresponds to the results obtained by Frank [16[ . 

For d — 3 and B — 2, equation ([5]) is simplified and an analytical solution can be found, resulting in 



u(z) ~ 



, e z2 / i (z 2 + 2) V^Ae A2 (erf(|)-erf(A) 



1 



2Ae 



-z 2 /4 



' A 3 e 



2z 



V7rAe A2 erfc(A) - 1 
: e z2 / 4 (z 2 +2)erfc(f) 



^ V^Ae A2 erfc(A) - 1 

Putting u(2\) = u s , we obtain 

l + 2A 2 (l- V^Ae A2 erfc(A) 



z 2 /4 
-n—lu- 



k = 



2A 2 -(l + 2A 2 )V^Ae A erfc(A) 



2(v^Ae A2 erfc(A) 



1 



2 A 2 (^y/TrX e A2 erfc(A) - 1 
For u s << Up, we write 

k = -2A 4 + V^A 3 (1 + 2A 2 ) e A2 erfc(A) + ( 1 - 2A 2 + 2V^A 3 e A2 erfc(A)) e + C(e 2 ). 



(20) 



(21) 



(22) 



General analytical expressions of u(z) and k, in terms of confluent hypergeometric functions, can also be found 
for even values of B as detailed in appendix El Furthermore, asymptotic forms of u(z) for large and small z can be 
calculated, see appendix [X] for analysis of z >> 1. Figure [3] compares k versus A for d = 2 and d = 3 in the absence 
of dislocation field (B = 0). 



IV. DISCUSSION 



The potential energy in equation ([3]) describes the elastic energy of the dislocation relaxed within the volume 
occupied by the second-phase precipitate [Ioj |. It was treated here as an external field affecting the diffusion- limited 
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FIG. 3: Growth coefficient A as a function of supersaturation parameter k at B — for a circular plate (d = 2) versus a sphere 
(d = 3). 



growth of second-phase precipitate. The interaction energy of impurities in a crystalline with dislocations depends 
on the specific model or configuration of a solute atom and a matrix which is used. Commonly, it is assumed that 
the solute acts as an elastic center of dilatation. It is a fictitious sphere of radius R' embedded concentrically in a 
spherical hole of radius R cut in the matrix. If the elastic constants of the solute and matrix are the same, the work 
done in inserting the atom in the presence of dislocation is w — pAv, where p is the hydrostatic pressure and Aw is 
the difference between the volume of the hole in the matrix and the sphere of the fictitious impurity. For a screw 
dislocation p = 0, while near an edge dislocation p = ^ 1 ^\ bG J) 1 " 6> f° r an impurity with polar coordinates (r, 9) with 
respect to the dislocation Oz, hence w oc Avsin9/r [Hi]. Using a nonlinear elastic theory [26j], a screw dislocation 
may also interact with the spherical impurity with the interaction energy w oc Av/r 2 . Moreover, accounting for 
the differences in the elastic constants of a solute and a matrix, the solute will relieve shear strain energy as well 
as dilatation energy, which will also interact with a screw dislocation with a potential w oc Av/r 2 [24]. Indeed, 
Friedel [24| has formulated that by introducing a dislocation into a solid solution of uniform concentration c , the 
interaction energy between the dislocation and solute atoms can be written as w ±2 wo(b/S) n f(9), where 6 is the 
distance between the two defects, wq the binding energy when d = &, and f(6) accounts for the angular dependence 
of the interaction along the dislocation. Also, n — 1 for size effects and n — 2 for effects due to differences in elastic 
constants. The discussed model for the interaction energy between solute atoms and dislocations has been used to 
study the precipitation process on dislocations by number of workers in the past [2^ . |28| and thoroughly reviewed in 
[TH . These studies concern primarily the overall phase transformation (precipitation of a new phase) rather than the 
growth of a new phase considered in our note. That is, they used different boundary conditions as compared to the 
ones used here. 

Let us now link the supersaturation parameter k to an experimental situation. For this purpose, the values of u s , 
i.e. the concentration at the interface between the second-phase and matrix should be known. The capillary effect 
leads to a relationship between u s and the equilibrium composition u eq (solubility line in a phase diagram) . To obtain 
this relationship, we consider an incoherent nucleation of second-phase on a dislocation a la Cahn [10(. A Burgers loop 
around the dislocation in the matrix material around the incoherent second-phase (circular plate) will have a closure 
mismatch equal to b. Following Cahn, on forming the incoherent plate of radius R, the total free energy change per 
unit length is 

g = -nR 2 Ag v + 27r 7 i? - A' ln(i?/r ), (23) 
where Ag v is the volume free energy of formation, 7 the interfacial energy and the last term is the dislocation energy, 
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A 1 = Gb 2 /An for screw dislocations, cf. equation ([3]). Setting dQ/dR = 0, yields 



R =2kS 1±VT ^)> (24) 

where a = 2A' Ag v /ir~f 2 . So, if a > 1, the nucleation is barrierless, i.e., the phase transition kinetics is only governed 
by growth kinetics, which is the subject of our investigation here. If, however, a < 1, there is an energy barrier and 
the local minimum of Q at R = Rq, which corresponds to the negative sign in equation (|24p , ensued by a maximum at 
R = R* corresponding to the positive sign in this equation. The local minimum corresponds to a subcritcal metastable 
particle of the second-phase surrounding the dislocation line, and it is similar to the Cottrell atmosphere of solute 
atoms in a segregation problem. When a — 0, corresponding to B = 0, the two phases are in equilibrium and the 
maximum in Q is infinite, as for homogeneous nucleation. 

For a dilute regular solution, Ag v = (ksT/Vp) h\(u s /u eq ), where V p is the atomic volume of the precipitate com- 
pound, u s is the concentration of the matrix at a curved particle/matrix interface and u eq that of a flat interface, 
which is in equilibrium with the solute concentration in the matrix. Equation (|24[) gives Ag v = j/R— A'/2ttR 2 . 
Hence, for a dilute regular solution, we write 

■c 

-RV R, 



U eq exp 



s^-s)]- < 25 » 



(26) 



where £ = /3V p j, (3 = 1/ksT and r\ = A 1 /2n^. Subsequently, the supersaturation parameter is expressed by 

, _ M eg exp[-£(1 - §)] - u m 
u p - M eg exp[-^(1 - §)] 

Taking the following typical values: 7 = 0.2 Jm~ 2 , G = 40 GPa, and b = 0.25 nm, then A' w 2.0 x 10~ 10 N and 
i] = 0.16 nm. Figure [5] depicts u s /u eq , from equation (f2"5| . as a function of scaled radius R/( for V p = 1.66 x 10 -29 
m 3 , rj = and 77 = 0.16 nm at T — 600 K. Equation ([23)) is analogous to the Gibbs-Thomson-Freundlich relationship 
[3| comprising a dislocation defect. 

Recalling now the values used for the interaction parameter B in the computations presented in the foregoing 
section, we note that for B = 2 and the above numerical values for G and b at T — 1000 K, we find I w 0.14 nm, 
which is close to the calculated value of rj. 

In Calm's model, the assumption that all the strain energy of the dislocation within the volume occupied by the 
nucleus can be relaxed to zero demands that the nucleus is incoherent. For a coherent nucleus forming on or in 
proximity of dislocations, this supposition is not true. Instead, it is necessary to calculate the elastic interaction 
energy between the nucleus and the matrix, which for an edge dislocation is in the form G& 2 /[47r(l — v)r\ for the 
energy density per unit length [2^]. In the same manner, to extend our calculations for growth of coherent precipitate, 
we must employ this kind of potential energy, i.e. the potential energy of the form <f>(r) = — Aln(r/ro) + C sin#/r, in 
the governing kinetic equation rather than relation ([5]). 



APPENDIX A: EVALUATION OF SOLUTIONS OF EQUATION © FOR d = 3 

For an ordinary second-order differential equation with a regular singularity, the Frobenius method can be used to 
obtain power series solution. On the other hand, when singularity is irregular, no convergent solution may be found; 
nevertheless, albeit divergent, the solution can be asymptotic. Let us write equation ([5]) for d = 3 in a generic form 

u" +p(z)u +q(z)u = 0, (Al) 

where primes denote differentiation with respect to z, p(z) = z/2 + (2 + B)/z and q(z) = B/z 2 . Since we have 
imposed the boundary condition u(oo) = u m , it is worthwhile to explore the behavior of the solution as z — > 00. 
But, first let us put equation (|A1|) in a more convenient form by setting u(z) = w(z)exp[— i J p(z)dz], which gives 
u" + (q — p'/2 + p A )u = 0. Here, without loss of generality, we consider 

u" + f(z)u = 0, (A2) 
B + 3 B(B-2) 
16 ~ 4^2 

Since f{z) is not 0(z~ 2 ) as z — > 00, then the point at infinity is an irregular singularity for u{z). We now look for 
solutions of (|A2J> by considering 



where /(*) = -- - -I L> (A3) 



u(*)~exp(53Vn(2)), (A4) 



71 = 
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FIG. 4: The size dependence of the concentration at the curved precipitate/matrix interface u 3 relative to that of the fiat 
interface u eq for a set of parameter values given in the text, cf. eq. (|25[) . 



TABLE I: Solutions to equation (IA5|I . 



Sequence Solution 1 Solution 2 



ipo 


z*/8 


-z?/8 


i'l 








1p2 


(5±2) In z 


-(S±i)lnz 


1p3 








1p4 


Bz~ 2 


~3(B + 2)z- 2 



where {il> n (z)}, n = 0, 1, . . . , is an asymptotic sequence as z — > oo. Substituting (I A4|) into equation (|A2|) 

< + V'i + • • • + (V-o + + • ■ ■ ) 2 + f(z) ~ 0, (A5) 

where we have tacitly assumed that ^J n {z) is (twice) differentiable and the resulting series are still asymptotic. 

Equation (|A5|) is used to determine the {i/v^z)}, n = 0, 1, . . . by successively applying the asymptotic limit z — > oo. 
The results for the first few terms are shown in table HI Hence, we write for z — > oo: 



B 



.(z) ~ A ie z / 8 z^ 1+B ^ l + ii + 0(z~ 4 ) 



w_ (z) 



A 2 e- Z 2 /8 z -(B/2+2) 



3(B 



+ G(z- 4 ) 



(A6) 
(A7) 



where A\ and A 2 are arbitrary constants. Note that the solution (|A6[) is divergent for large z, whereas (|A7|) is 
convergent and thus is physically admissible. Considering w(oo) = u m , we write 



U-(-z) ~ u m + A 2 e 



r-B/2-2 



3(5 + 2) +C(z _ 4) 



as z — > oo. 



(A8) 



Let us now evaluate the general solution to equation ([5]) for d = 3 as expressed by equation (J9j) . We apply the flux 
conservation relation ^ to obtain C\, and then substitute C\ in equation @ to write 



u(z) =K 1 {z,B)C 2 +K 2 (z,B)q, 



(A9) 
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TABLE II: Special cases of iF 1 (a;b; z). 

1 F 1 (-l;±;-z*) = l + 2z* 

i-Fi(-|;|;-« 2 ) = e- z2 + v /5Fzerf(z) 

xFiC-i; |; -z 2 ) = + ^(1 + 2 2 2 )erf (^z" 1 

iF 1 (|;|;-» 2 ) = ^erf(«)z- 1 

f;- 2 2 ) = f (W^ + y^(2 2 2 ~l)erf(z)) Z - 3 



where 



Ki(z,B) 



-) ,/•)( - 



2 



+ 



R\-B Z? f 2-B . 5-B . 

B-3 



K 2 (z,B), 



K 2 (z,B) 



K 2 {z,0) 



4( J B + l)A 3 1 F 1 (-i;l±^; 



(B* - l)iF!(-|; ±±^; -^ 2 ) + 2A 2 iF 1 (|; -A*) " 
r2e- z2 /4 



2A 3 e A 



%erf(z/2) 



(A10) 
(All) 

(A12) 



Here, %Fi (a; 6; z) is the confluent hypergeometric function. If a < 0, and either b > or b < a, this function can 
be expressed as a polynomial with finite number of terms. If, however, b — or a negative integer, then 1 i 7 ' 1 (a;6;2;) 
itself is infinite. Thus, relations (|A11[) - (|A12[) become singular for B = 1, 3, 5, 7, . . . , making the solutions meaningless. 
Some useful relations for computations are listed in table [Til Additional relations and properties for iF 1 (a;b] z) can 
be found in [25j . 

Next, we utilize the remote boundary condition u(oo) = u m to determine C 2 , then we formally write 



u(z) 



Kx{oo,B)' 



(K 2 {z,B) 



^1(00,5) 



K 2 (oo,B))q 



(A13) 



In computations of Ki(oo, B) prudence must be exercised, i.e., first evaluate this quantity for a given value of B, then 
take the limit z — > 00. Note also that Ki(z, 0) = 1 Vz and ^2(00, 0) = 2y / 7rA 3 e A . 

Furthermore, we may calculate a relation for the supersaturation parameter, k = (u s — u m )/{u p — u s ), defined in 
the main text by using the condition u(2A) = u s on equation (|A9|) . which gives 



For dilute alloys, u s << u p ; so with e = u s /u p , we write 



q 2V ' ' if 1(00, £?) V ; 



(A14) 



fc 



K X {2\,B)) Kx{2\,B) u ' ; v ; v y 



(A15) 
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